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Abstract— Todo sistema eléctrico de potencia puede 
considerarse como formado por un conjunto de 
cuadripolos, que son los elementos básicos del sistema: 
en efecto, un transformador, una línea, un cable, una 
impedancia cualquiera, una carga pasiva, etc., pueden 
todos representarse mediante un cuadripolo pasivo y por lo 
tanto el cálculo de la red requiere el conocimiento general 
del funcionamiento del cuadripolo y de sus algoritmos de 
cálculo. La red es en definitiva una combinación de 
cuadripolos. 

El objetivo de este trabajo es el de presentar un método 
general de cálculo para la transferencia de potencia a 
través de un cuadripolo. Este método permite calcular, a 
partir de 3 datos eléctricos, todas las demás magnitudes 
que intervienen en dicha transferencia. 

Para la realización de ese objetivo, se presenta en este 
trabajo un programa de cálculo que hemos denominado 
"RESCUAD" y que se basa en el método de Newton - 
Raphson para la resolución de ecuaciones no lineales. 


Keywords— Cuadripolo, régimen alterno sinusoidal, 
transferencia de potencia, método de Newton-Raphson. 


I. INTRODUCCIÓN 


Es bien conocido que en un sistema eléctrico de 
potencia funcionando en régimen alterno sinusoidal, 
el modelo básico que aparece en toda porción pasiva 
del circuito representada por fase es el llamado 
“cuadripolo”, que tiene 2 bornes de entrada y 2 de 
salida: 





Fig. 1: Representación del cuadripolo. 
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Las convenciones usuales para las tensiones U, y 


U de entrada y salida respectivamente, así como 
para las corrientes /, e Z, son las indicadas en la 


figura. 

En todo este trabajo, las magnitudes señaladas 
con letras suprarayadas son fasores complejos o 
números complejos; para una magnitud cualquiera, si 
le suprimimos la raya superior, estaremos indicando 
el módulo de esa magnitud o de ese complejo; así por 
ejemplo: 


V fasor complejo; V= mód(V ) = V| 


Para el cuadripolo CUA representado en la 
figura, el borne asociado al signo + será la fase, 
mientras que el asociado al - será el neutro del 
sistema. Las tensiones que hemos indicado son 
entonces tensiones “estrelladas”, o sea tensiones de 
fase respecto al neutro, y las corrientes serán 
corrientes de línea. 


Todo elemento pasivo en un circuito, sea una 
impedancia, sea una línea, sea una combinación de 
impedancias, puede reducirse a un cuadripolo. 
También cuadripolos en paralelo (caso de líneas en 
paralelo) o en cascada pueden reemplazarse por un 
único cuadripolo equivalente, empleando las 
conocidas fórmulas de álgebra de cuadripolos 
deducidas de operaciones matriciales. 


Las representaciones más usuales del 
cuadripolo de potencia en la teoría de los sistemas 
eléctricos de potencia son el modelo en “TI* y el 
modelo en “T”, o también más frecuentemente el 
modelo de matriz cuadrada de constantes generales: 


A 
T — > 


siendo también conocidas las fórmulas de pasaje de 
un modelo a otro. 

En este trabajo , nuestro cuadripolo CUA estará 
dado por su matriz de constantes generales: 


A=Ae”” (adimensionada) 
B= Be” (Q) 
C=Ce* (Q) 


D= De”? (adimensionada) 


En el cuadripolo pasivo de potencia, estas constantes 
no son independientes pues están ligadas por la 


relación AD -— BC =1, como lo muestra la Teoría 
de Circuitos. Será entonces suficiente conocer 3 de 
esas constantes para que el cuadripolo esté 
determinado. Más aún, si el cuadripolo es simétrico 
(por ejemplo el caso de una línea de transmisión) se 


tiene D = A , de modo que en ese caso alcanzará el 
conocimiento de 2 de esas constantes. 

Las constantes generales del cuadripolo 
intervienen en el problema eléctrico de acuerdo a la 
ecuación matricial: 

Uı A B||U2 


n| |c Dp]||r 
o sea: 

Ui =AU2+BI, 

Iı =CU2+DI> 


Observemos que estas 2 ecuaciones en complejos 
equivalen a 4 ecuaciones con magnitudes escalares, 
en las que intervienen 8 variables (módulos y 
argumentos de 2 tensiones y 2 corrientes). Pero 
haciendo notar que el origen de fases puede ser 
libremente elegido, podemos fijar en cero uno 
cualquiera de los argumentos y por lo tanto quedan 
sólo 7 variables. 

4 ecuaciones con 7 variables nos dejan 3 grados 
de libertad. Para que un problema de 
transferencia de carga a través de un cuadripolo esté 
determinado, será necesario y suficiente conocer 3 
magnitudes eléctricas escalares, que constituirán 
siempre el punto de partida del problema. 





Il. PLANTEO GENERAL 


Trabajaremos con las magnitudes usuales en los 
cálculos de sistemas eléctricos de potencia, o sea: 
- tensión 
- potencia activa 
- potencia reactiva - 
factor de potencia 


Teniendo en cuenta los dos lados del cuadripolo, 
tenemos pues 8 magnitudes en juego, para las cuales 
tomaremos las siguientes convenciones de sentido: 


+ + 
— — 
B ¡Papo á 
2 


Fig. 2: Convenciones de sentido en las potencias y 
tensiones del cuadripolo. 


Las potencias aparentes se han tomado entrante para 
la entrada y saliente para la salida; esas potencias 
aparentes son: 


Sı =P, + jQ, 
S2 =P,+j0, 


Los factores de potencia a la entrada y a la salida se 
evalúan a través de los argumentos de esas potencias 
aparentes, o sea: 


ọ = arg(S1) 
p, =arg(S2) 


Las 8 magnitudes que intervendrán en nuestro 
problema son pues: U;, Pi, Qi, (1, U2, P2, Q2, Q» . 
En muchos problemas de transferencia de potencia, 
resulta útil conocer también el ángulo entre las 


tensiones terminales U, y U2, por lo cual 
pondremos en juego también el adelanto O de U, 


respecto a U2 , quedándonos en definitiva un juego 
de 9 magnitudes escalares: U1, P1, Q1, (1, U2, P2, Q», 
02, 6. 

Dado que trabajamos por fase, las potencias P y Q que 
consideramos deben tomarse también por fase (o sea 
la tercera parte de las potencias trifásicas que 
intervienen). 


Vistos los 3 grados de libertad señalados 
anteriormente, tendremos que escribir 6 ecuaciones 
que relacionan a esas 9 variables. 


HI. PLANTEO DE LAS ECUACIONES 


A. Introducción 

Dejando de lado moméntaneamente el ángulo 
0, tendremos que escribir 5 ecuaciones, o sea 2 
ecuaciones en complejos más la condición 


mencionada de origen de fases en una tensión o una 
corriente. 


Tomando por ejemplo U¡ como origen de 


fases U¡ =U 1» las dos ecuaciones serán: 
o =£(0,5) 0) 
S2 = f,(U,,81) (2) 
ya que resulta lógico pensar en calcular la condición 
eléctrica de salida en función de la condición eléctrica 
de entrada. 


B. Obtención de la ecuación (1) 


La inversión de la ecuación matricial básica nos 


BEJE 


habiendo tenido en cuenta que AD-BC=1. 


Tenemos pues: 


U2 =DU,-BI, 
Reemplazando h= S (el símbolo “^” significa 
conjugación): 
U2 =DU,- pá (1) 
U 


1 


C. Obtención de la ecuación (2) 


Tenemos: 
S2=U2Î, 


Siendo 12 = -CU i +All, , tendremos 


Î, =-CU,+AÍ, =-CU, +AD , de donde: 
1 


A 


S2 = (DU, BÍO -ĈU + 2 
U, U, 


Operando, se llega fácilmente a: 





7 e A ÂB 
S2=ADS¡+BCS, -ĈDU ? ngas (2) 


1 
D. Forma escalar de las ecuaciones (1) y (2) 


Igualando los cuadrados de los módulos de ambos 
miembros de la ecuación (1), se llega a: 


P? + 2 
U; =D’°U? +B? 22 2 - o 
U: a’) 


2BDÍP, cos(B —8)+0, sen(ß -ô)] 


Igualando componentes reales e imaginarias de 
ambos miembros de la ecuación (2), se llega a: 


P, = ADÍP, cos(æ-8)+ Q, sen(& -8)]+ 
BC[P, cos(B—y)+0, sen(B=y)]- 

2 2 
e Hi cos(f -&) 


1 


- CDU cos(ô —y)- AB 





(2r) 


Q, = ADÍP, sen(&- 8) +Q, cos(æ -8)]+ 
BCÍ[P, sen(PB -8) +Q, cos(8 -8)]- 

2 2 
E o sen(P —02) 





CDU | sen(ô —y)- AB 
(21) 
E. Ecuaciones que hacen intervenir Q, y Q> 


Estas son, obviamente: 


P¡senQ, = Q¡cosp; 6) 
P>senQ, = Q,cosq, (4) 


F. Ecuación para el adelanto O de U, respecto a 


U, 


Habiendo tomado U¡ = U ı como origen de fases, 


tendremos O = -argU2. La ecuación (1) puede 
entonces escribirse: 


o Poj 
Ue” = DeP*G, — Be? LB 41 jQ, , 
U, 
o sea: 


U,Uze* = DU/%e”-B(P,-¡Qye? 


Igualando componentes reales e imaginarias de 
ambos miembros, se llega fácilmente a: 


U¡U,cos0 = DU¡?cosó - B(P,cosB+Q¡senf) (5) 
U¡U,sen0 = -DU/*send + B(P¡senf-Q¡cosB) (6) 


Para determinar 0, usaremos sólo la ecuación (5), 
que nos da 2 valores opuestos de O a partir de cosb; 
el valor correcto a introducir en la solución se elige 


de acuerdo al signo de sen dado por la ecuación (6). 


G. Conjunto de ecuaciones a resolver 


Hacemos un listado de las ecuaciones halladas, 
reformulándolas para tener un segundo miembro 
nulo: 


UŻU} -D°’Uf -B° (P? +0? )+ 
2BDÍ?, cos(P-8)+0, sen(P—8)|=0 
(ec.1) 


P, — ADÍP, cos(& -8)+ Q, sen(œ - 8)]- 

BCÍP, cos(B—y)+0, sen(B—y)]+ 

P? +O 
1 


+CDU? cos(ô — y) + AB cos( PB —0) =0 


2 


(ec.2) 


Q, — ADÍP, sen(a—8)+0, cosía -9)]- 
BC|P, sen(B=8)+0, cos(P -8)]+ 
E +0, 


2 
1 


+ CDU sen(9 —y)+ AB sen(B=0)=0 
(ec.3) 


P¡senp, - Q¡cosp, = 0 
(ec.4) 


P>seno, - Q2cosp, = 0 
(ec.5) 


U¡U,cos8 - DU, cos - B(P,cosB+Q¡senf) = 0 
(ec.6) 


Como puede verse, tenemos las 9 variables: 


Us, Pi, Qi, Q1, U2, P2, Q2, Q», O 


con 6 ecuaciones, o sea 3 grados de libertad como 
señalado anteriormente. 

El sistema obtenido se resuelve por iteración 
utilizando el método de Newton-Raphson. Se elige 
además la determinación correcta de 0 verificando 
que: 


U¡U,sen0 +DU“send - B(P¡senf-Q¡cosB) = 0 


Notas : 
1. El problema planteado se refiere a un sistema 
trifásico equilibrado, en que los 2 polos del 
cuadripolo en cada terminal son una fase y el neutro 
del sistema. Por lo tanto, las tensiones son tensiones 
estrelladas y las potencias son potencias por fase. Una 
simple observación aritmética muestra que pueden 
utilizarse las mismas fórmulas con tensiones 
compuestas y potencias totales. El programa 
propuesto presenta pues la ventaja de poder trabajar 
directamente con las magnitudes manejadas 
usualmente en los problemas de sistemas trifásicos 
equilibrados. 
2. La eliminación gaussiana incluida en el 
método de Newton-Raphson es válida siempre que el 
sistema de ecuaciones cuya solución da el paso 
siguiente para la iteración sea determinado: Si el 
sistema es indeterminado, el programa acusará 
“error”. Un ejemplo trivial en que el sistema sería 
indeterminado es el caso en que los 3 datos sean los 
siguientes: 

U; Q:=0 @=0 , 
que correspondería a conectar en la salida una 
resistencia pura R. Es fácil comprender que el sistema 
en este caso es indeterminado pues el valor de R no 
aparece en los datos y si el método nos diera una 
solución, llegaríamos al absurdo de que ¡¡el 
funcionamiento es independiente de R!! El programa 
realizado envía un mensaje de error invitando al 
usuario a cambiar los datos iniciales (ver Apéndice). 


IV. PROGRAMA RESCUAD 


Este programa resuelve el sistema de 
ecuaciones por iteración empleando, como ya 
mencionado, el método de Newton-Raphson. 

La entrada y salida de datos se efectúa en forma 
directa a través de la pantalla del PC. Se optó por esta 
modalidad en vez del ingreso de datos en un archivo 
por considerar que la misma permite una mejor 
interacción del usuario con el programa 
implementado. 

A continuación se muestran las pantallas 
típicas de entrada y salida para un caso particular: 


Esa línea está alimentada por una red que le impone 
una tensión de 380 kV a la entrada y le entrega una 
potencia de 200 MW; a la salida, la línea alimenta 


380 kV una carga sélfica de cosp=0,8. Se desea conocer todas 
> L > las magnitudes eléctricas en la entrada y en la salida 

Red. Carga de la línea. 

200 MW sélfica Para el programa RESCUAD, tenemos el 


cosq=0.8 cuadripolo equivalente a la línea, simétrico y de 


constantes A,B conocidas, con los siguientes datos e 


incógnitas: 
e 
Se : 5 j p P, = 200 MW 
equivalente: 
, = Arccos0,8 = 36,87" 
0,9271 e" Incógnitas :Q,, Q1, U2 , P2 , Q2, 8 


wia] 
| 


NON Mostraremos lo que aparece en pantalla: 
96,23 e*2% (Q) 


PROGRAMA DE CALCULO DE CUADRIPOLOS IIE 
- Facultad de Ingeniería - 
INGRESO DE LOS DATOS DEL PROBLEMA 


Ingrese la tolerancia deseada: 0.01 
Ingrese el número máximo de iteraciones admitidas: 30 


Constantes del cuadripolo: 


Indique si el cuadripolo es simétrico (S/N) S 


Ingrese el valor de A:0.9271 
Ingrese el valor de alfa(grados):0.791 


Ingrese el valor de B(ohm):96.23 
Ingrese el valor de beta(grados):82.66 
Variables del problema 
Puede elegir dos sistemas de unidades: 
1) [UJ=V, [P]=W, [QI=VAR 
2) [U]=kV, [P]=MW, [Q]=MVAR 


U1: Es dato? (S/N) S 
Valor: 380 


P1: Es dato? (S/N) S 
Valor: 200 








Q1: Es dato? (S/N) N 
Valor: 150 


FI1: Es dato? (S/N) N 





Valor (en grados): 36.87 


U2: Es dato? (S/N) N 
Valor: 380 


P2: Es dato? (S/N) N 
Valor: 200 


Q2: Es dato? (S/N) 
Valor: 150 


FI: Es dato? (S/N) S 
Valor (en grados): 36.87 


theta: Es dato? (S/N) N 
Valor(en grados): 10 


Los datos ingresados fueron: 
Constantes del cuadripolo: 


A: 9.2710000000E-01 alfa: 7.9100000000E-01 
B: 9.6230000000E+01 beta: 8.2660000000E+01 
C: 1.4839298060E-03 gamma: 8.7774671922E+01 
D: 9.2710000000E-01 delta: 7.9100000000E-01 


Datos y valores iniciales de las incognitas: 


U1: 3.8000000000E+02 S 
P1: 2.0000000000E+02 S 
Q1: 1.5000000000E+02 N 
FI1: 3.6870000000E+01 N 
U2: 3.8000000000E+02 N 
P2: 2.0000000000E+02 N 
Q2: 1.5000000000E+02 N 
FI2: 3.6870000000E+01 S 
theta: 1.0000000000E+01 N 





Nota: Si al correr el programa, en alguna de las iteraciones el sistema 
de ecuaciones resulta incompatible o indeterminado, verifiquesusdatos 
o modifique ligeramente los valores iniciales de sus incógnitas. 


RESULTADOS: 


Q1: -4.2105608180E+01 
FI1: 3.4811124731E+02 
U2: 3.5942007149E+02 
P2: 1.8678109527E+02 
Q2: 1.4008634281E+02 
theta: 7.4706337110E+00 
niter: 10 
converge: TRUE 





Con 10 iteraciones, el programa nos brinda 
todas las magnitudes eléctricas en el 
funcionamiento planteado: 


Q; = -42,11 MVAR 
Ọı = 348,1° (o sea -11,9°) 


U? = 359,4 kV 

P2 = 186,8 MW 
Q2 = 140,1 MVAR 
0 = 7,47° 


La red recibe potencia reactiva (la línea es muy 
capacitiva), las pérdidas en la línea son 200 - 186,8 
= 13,2 MW, la caída de tensión en la línea es del 


5,42% y el adelanto de Us respecto a U? es de 
7,47". 


V. OBSERVACIÓN SOBRE EL INGRESO DE 
DATOS EN EL PROGRAMA RESCUAD 


Para acelerar la convergencia del cálculo, se 
recomienda tomar como valores iniciales de las 
incógnitas valores próximos a los de las variables 
homólogas conocidas: por ejemplo si se conoce P4 
y no se conoce P,, tomar como valor inicial de P, el 
dato P4; si no se conocen las 2 variables homólogas 
Qı, Qz, tomar valores compatibles con el orden de 
magnitud fijado por ejemplo por P,, Q1. Debe 
tenerse en cuenta que en un sistema eléctrico de 
potencia, el cuadripolo es simplemente un 
elemento de tránsito de carga y que no modifica 
substancialmente las variables eléctricas de la 
entrada, produciendo pues en la salida valores 
cercanos a los de la entrada, fundamentalmente en 
lo que se refiere a tensiones y potencias activas. 
Análogamente, para el ángulo O conviene 


comenzar la iteración suponiendo U' yU, en fase, 


o sea por ejemplo comenzar con 9 =0. 
En el ejemplo mostrado se tomó O = 10° como 
valor inicial. 


VI. APENDICE — MENSAJE DE ERROR 


El programa RESCUAD no conduce a solución 
alguna cuando se produce una de estas 2 
situaciones: 


A) Los datos del problema son incompatibles o 
redundantes. 

B) Los valores iniciales elegidos para las 
incógnitas producen en alguna de las etapas de la 
iteración un sistema de ecuaciones indeterminado o 


incompatible (o sea un sistema con determinante 
principal nulo). 

Antes de correr el programa, la pantalla exhibe un 
mensaje al ususario: 


"Si al correr el programa, en alguna 
de las iteraciones el sistema de 
ecuaciones resulta incompatible o 


indeterminado, verifique sus datos o 
modifique ligeramente los valores 
iniciales de las incógnitas" 





Luego, al correrlo, si se produce una de las dos 
situaciones mencionadas, el programa exhibe el 
mensaje: 


"HORROR! Sistema de ecuaciones 


incompatible o indeterminado" 





Veamos algunos ejemplos triviales en que se 
producen estas situaciones: 


Situación A : 


A1) Ejemplos de datos incompatibles 





Ejemplo A.1.1 
Nos damos los 3 datos P2, Q2, @2, con tgQ, + £ ; 
2 


Es evidente que este caso conduce a una situación 
de incompatibilidad. 


Ejemplo A.1.2 


Para un cuadripolo en que las constantes generales 
son todas imaginarias puras, o sea que no hay 
resistencias, no existe pérdida de potencia activa en 
el cuadripolo. Si entre los 3 datos figuran P, y P2 
con P; + P», es evidente que la situación resultará 
incompatible. 


A2) Ejemplos de datos redundantes 





Ejemplo A.2.1 
0, 


Nos damos los 3 datos P», Q2, (2, con tg, = po 
2 


Es evidente que uno de los 3 datos es redundante y 
que el sistema resultará indeterminado pues en 
realidad sólo hemos impuesto 2 condiciones 
eléctricas. 


Ejemplo A.2.2 


En el caso de una resistencia pura R conectada a la 
salida del cuadripolo, si nos damos por ejemplo U», 
Q- = 0, @ = 0, el sistema resultará indeterminado. 
Es previsible la indeterminación, pues si el 
programa suministrara una solución, ésta sería 
independiente de R, lo cual es absurdo. La 
indeterminación surge del hecho de la redundancia 
de las informaciones Q; = 0, q, = 0. En realidad 
nos estamos dando sólo 2 informaciones. 


Situación B 


Ejemplo de mala elección de valores iniciales 





Las ecuaciones planteadas muestran que si A = 1 
y P = 90°, tomar el valor inicial 9 = 0 conduce a 
un sistema indeterminado. Conviene partir de un 
valor O + 0 para poder llegar a la solución. 
Observemos que la situación de determinante nulo 
puede darse en la primera etapa de la iteración o en 
cualquiera de las etapas siguientes. El mensaje 
emitido por el programa no especifica en cual de 
las etapas se produce la anulación del deteminante, 
pero, de cualquier modo, si los datos del problema 
son adecuados, debe modificarse el punto de 
partida de la iteración cambiando (aunque sea 
ligeramente) los valores iniciales de las incógnitas 
o de alguna de ellas. 
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